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We study the strong and electromagnetic decay properties of scalar mesons above 1 GeV within 
a chiral approach. The scalar-isoscalar states are treated as mixed states of quarkonia and glueball 
configurations. A fit to the experimental mass and decay rates listed by the Particle Data group 
is performed to extract phenomenological constraints on the nature of the scalar resonances and to 
the issue of the glueball decays. A comparison to other experimental results and to other theoretical 
' approaches in the scalar meson sector is discussed. 
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I. INTRODUCTION 



> 

The unique interpretation of scalar mesons constitutes an unsolved problem of hadronic QCD. Below the mass scale 
CN * of 2 GeV various scalar states Q] are encountered: the isoscalar resonances a = /o(400 — 1200), /o(980), /o(1370), 
/o(150Q) and / (1710), the isovectors a (980) and a (1450) and the isodoublets K*(800) and K*(U30). The existence 
of the .Kg (800) is still not well established and omitted from the summary tables of |lj. From a theoretical point of 



view one expects the scalar quark-antiquark ground state nonet ++ , a scalar-isoscalar glueball, which lattice QCD 
£p-t predicts to be the lightest gluonic meson with a mass between 1.4-1.8 GeV |2j, and possibly other exotic states (non 
qq states), like e.g. four quark states or mesonic molecules Various interpretations of and assignments for the 
■ physical scalar resonances in terms of the expected theoretical states have been proposed (see, for instance, the review 
. ^ ' papers [J Q Q and references therein) . 

In this work we follow the original assignment of Ref. where in a minimal scenario the bare quarkonia states 

h ; 

N = ^l/2(uu + dd) = nn, S = ss and the bare scalar glueball G mix, resulting in the three scalar-isoscalar resonances 
/o(1370), /o(1500) and /o(1710V. Such a mixing scheme has been previously investigated by many authors, as for 
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example in the lattice study of [jj or within the model approaches of 

The mesons cio(1450) and A"g(1430) are considered as the 1=1 and 1=1/2 quarkonia J PC = ++ states. In 
this way the low-lying scalar quarkonia nonet is located in the energy range of 1-2 GeV, where the other p-wave 
nonets of tensor (2 ++ ) and pseudovector mesons (1 ++ ) Q are also situated. Masses of selected scalar quarkonia 
states were also estimated on the Lattice; in Ref. |14| the 1=1 scalar quarkonium state is predicted to have a mass 
of M ao = 1.51 ± 0.19. This result favours the interpretation of the state arj(1450) (and not ao(980)) as the isovector 
ground-state scalar quarkonium. However, previous lattice studies (as in |l5| : see also |l llllul and Refs. therein) find 
different results. Scalar resonances below 1 GeV can possibly be interpreted as four quark states or mesonic molecules. 
In the hypothesis of Jaffe's four quark states was studied on the Lattice, where the masses of four-quark states 
are found to be lighter than 1 GeV. Recent attempts in the context of chiral perturbation theory to describe the 
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scalar states below 1 GeV as "dynamically generated" resonances, i.e. states which do not survive in the large N c 
limit Q ! have also been performed. In this scheme the nature of the scalars below 1 GeV can also be related 
to four quark configurations. This can be viewed as a further indication that scalar quarkonia states are located at 
masses above 1 GeV. 

In this work, starting from an effective chiral Lagrangian derived in Chiral Perturbation Theory (ChPT) [3, 
I21I l^ . we perform a tree- level analysis of the strong and electromagnetic decays of scalar mesons settled in the 
energy range between 1 and 2 GeV. The scalar glueball is introduced as an extra-flavor singlet composite field with 
independent couplings to pseudoscalar mesons (and to photons, although suppressed). Although a chiral approach 
cannot be rigorously justified at this energy scale, since loop corrections could be large, we intend to use this framework 
as a phenomenological tool to extract possible glueball-quarkonia mixing scenarios from the observed decays. 

The scalar glueball G mixes with the scalar quarkonia fields N = ^J\j1{uu + dd) = nn and S = ss in accord with 
flavor blindness. We also consider in this scheme a possible direct mixing of the quarkonia fields N and S. The origin 
of such mixing can be driven by instantons [2^ . The presence of a (even small) flavor mixing in the scalar- isoscalar 
sector can sensibly affect the phenomenology. 

In the presented approach the glueball decay into two pseudoscalar mesons is occurring by two mechanisms: 
a) through mixing, that is the glueball G acquires a quarkonium component, which subsequently decays into two 
pseudoscalars^; b) direct deca y o f the glueball component G into two pseudoscalars without an intermediate scalar 
quarkonium P, Q Q- In jr3. Ilf| this direct decay is argued to be suppressed as based on arguments of the strong 
coupling expansion, while in the phenomenological fit of Q] it dominates. Here we also address the question, if the 
direct decay is needed to explain the decay phenomenology of the scalar-isoscalar resonances /o(1370), /o(1500) and 

/o(mo). 

For consistency we also consider the strong decays of the isovector and isodoublet scalar states as well, trying to 
highlight the difficulties and the open issues, also comparing with previous works on this subject. Following the idea 
of |20j we analyze deviations from the large N c limit in the framework of the proposed scalar-quarkonia assignment 
(which differs from 20] ), which turn out to be relatively small in the phenomenology. 

The paper is organized as follows: in Section [H] we discuss the effective Lagrangian related to the scalar quarkonia- 
glueball mass spectrum and the strong and electromagnetic decays. In Section lLTTl we determine a phenomenological fit 
to the experimental data listed in yj by first neglecting the direct decay of the glueball component. In Section llVl we 
also allow for a direct glueball decay studying its influence on the results. There we use the lattice data of [2^, where 
an approximate calculation of the glueball decays into two pseudoscalar mesons has been performed, to constrain our 
analysis. Finally, in Section we summarize our results and draw conclusions. 

II. THE MODEL 
A. The Lagrangian 

The strong and electromagnetic decays of scalar mesons are based on an effective chiral Lagrangian C c g as derived 
in Chiral Perturbation Theory (ChPT) j^l^Q- The Lagrangian involves the nonets of pseudoscalar and of scalar 
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'■a^' 5 = 75^^., (1) 

V 1=0 V 2 — 

the electromagnetic field and, in addition, a new degree of freedom, the bare glueball field G, which is treated as a 
flavor-blind mesonic field. The lowest order effective Lagrangian C e g in the large N c limit including 1 /N c corrections 
reads 



.off = — (D^UD^ + X +) + ^(D^SD"S-MlS 2 ) + -(d^Gd fi G-M^G 2 ) 



+ c s d (S u„ «" ) + c s m (S X+ ) + -jf ( U/1 u") + G ( X+ ) 

+ cl (S F+ F+^)+^=G( F+ F+ + £? ix + £f lix . (2) 

Here the symbol (...) denotes the trace over flavor matrices. The constants c d , c s m , <A, c^, c\ and cf define the coupling 
of scalar fields and of the bare glueball to pseudoscalar mesons and photons, respectively. We use the standard notation 
for the basic blocks of the ChPT Lagrangian W'- U = u 2 = cxp(iVV2/F) is the chiral field collecting pseudoscalar 
fields in the exponential parametrization, denotes the chiral and gauge-invariant derivative, u M = iu' D^Uw is 
the chiral field, x± — v^x u ^ ± U X^ U -, X — 2-B(s + ip), s = M. + . . . and F& v = u'F^Qu + uF^Qu^ , where F^ v is 
the stress tensor of the electromagnetic field. The charge and the mass matrix of current quarks are denoted by Q — 
ediag{2/3, — 1/3, — 1/3} and M. = diagjm, to, m s } (we restrict to the isospin symmetry limit with m u — m d = to); 
B is the quark vacuum condensate parameter and F the pion decay constant. 

The masses of the octet pseudoscalar mesons in the leading order of the chiral expansion (first term of the La- 
grangian) are given by 

Ml = 2mB , M 2 K = (m + m s )B , M 2 8 = ~ (m + 2m s )B . (3) 

o 

The contribution to the mass of rj° in leading order is 

M 2 „ = ^(2m + m s )B , (4) 

i.e. rp is a Goldstone boson in the large iV c and in the chiral limits. 

Following [2^] , we encode in £„ lix an extra-contribution to the mass of rf (due to the axial anomaly) and the if -if 
mixing term: 

{f, )' - zprPrf , (5) 
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(the parameters "fp and zp are in turn related to the parameters M m and d m of |22j). The physical diagonal states 
77 and 77' are given by 

77° = rf cosOp — r\ sin#p , 77 s = rf sin 6>p + 77 cosOp , (6) 



where Op is the pseudoscalar mixing angle. We follow the standard procedure |2(J, l22t l25l |26| and diagonalize the 
corresponding rf-rf mass matrix to obtain the masses of 77 and 77'. By using M w — 139.57 MeV, Mk — 493.677 MeV 
(the physical charged pion and kaon masses), M n = 547.75 MeV and M v i — 957.78 MeV the mixing angle is 



determined as 9p = —9.95°, which corresponds to the tree-level result (see details in Ref. j^)- Correspondingly 
one finds M„o = 948.10 MeV and zp = —0.105 GeV 2 . Higher order corrections in ChPT cause a doubling of the 
absolute value of the pseudoscalar mixing angle j^); m our wol 'k we restrict to the tree- level evaluation, we therefore 
consistently use the corresponding tree-level result of 9p = —9.95°. In the present approach we do not include the 
neutral pion when considering mixing in the pseudoscalar sector, because we work in the isospin limit. This mixing is 
small, and can be safely neglected when studying the decay of scalar resonances into two pseudoscalars. Similarly, for 
all pseudoscalar mesons we use the unified leptonic decay constant F, which is identified with the pion decay constant 
F = F-n = 92.4 MeV. A more accurate analysis including higher orders should use the individual couplings of the 
pseudoscalar mesons (for a detailed discussion see Refs. |27|). 

B. Scalar quarkonia - glueball mixing 

In this subsection we discuss the glueball-quarkonia mixing. For this reason we restrict to the following part of the 
effective Lagrangian ©: 

C s = \ (D,SD»S - M 2 S S 2 ) + 1(8^0 - M 2 G 2 ) + (7) 

where, besides the scalar quarkonia nonet, also the scalar glueball field G has been introduced. In the large N c limit 
all the states of the quarkonia nonet have the same nonet mass Ms and the glueball is decoupled from the quarkonia 
sector. Deviations from this limit are encoded in £f lix , where the glueball-quarkonia mixing is introduced and where 
the degeneracy of the nonet states is lifted. 

For what concerns the explicit mass term and the next-to- leading order 1/N C terms in the quarkonia sector we 
follow 20]. Including additionally a possible breaking of the Gell-Mann-Okubo (GMO) mass relation and the glueball- 
quarkonia mixing under the hypothesis of flavor blindness, we have: 

M 2 M 2 

Aii* = e S m (S 2 X+ ) + k^So (x+> - 7s ~^S 2 - js 8 ^S% - VsfGSo . (8) 

The parameter e^, describes the strength of flavor symmetry breaking derived from the non-zero values of the current 
quark masses, the parameters and 7s describe the order 1/N C terms. The parameter 7s 8 and the related term 
in the Lagrangian is not derived in the iV c expansion (it is in fact absent in |20J), but it describes violations of the 
GMO mass formula (indeed a result from higher orders in the chiral expansion). Finally, the parameter / is the 
glueball-quarkonia mixing strength. Note that G, being a flavor singlet, couples only to the flavor singlet quarkonium 
state So in the flavor-blind mixing limit. The glueball and the quarkonia sector decouple in the large N c limit |3C)|: a 
non-zero mixing, as described by a non-vanishing parameter /, takes into account a possible deviation from the large 
N c limit, as the parameters k^ and 7s in the quarkonia sector. 

The terms contained in £^ ix lead to mass shifts of the nonet masses, while also introducing mixing both among 
quarkonia states and in the glueball-quarkonia sector. With the use of Eqs. J7J and JS| the explicit expression of the 
Lagrangian C s reads: 

C s = ~a (D + M 2 o )a K* * (□ + M£,) K* l -G{U + M 2 )G 

~ lM 2 8 S 2 ~hl 2 o S 2 ~z s SoSs-V3fGS (9) 



where ao is the isovector collecting and a[j fields; Kq and Kq are the doublets of {K^ + , ^o°) an< ^ C^o > -^0°) 
mesons, respectively (see Eji an d where the corresponding masses are given by (see Eqs. Q and JEJ)): 



By inverting we get: 



M 2 


= Ml-Ae s m Ml, 




M 2 K * 


= M%-Ae s m M 2 K , 






= Mj(l+ 7s .)-|e*(4M£ 


-M*). 




= M|(l+ 7So )-ie^(2M2 


+ M*). 
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(10) 



(11) 

We can immediately deduce Afs and for the considered assignment (which differs from j2(lj ) by using the experi- 
mental masses M ao = M ao ( 1450 ) = 1.474 ± 0.019 GeV and Mr* = Mk*(uw) = 1.416 ±0.006 GeV Q: 

Ms = 1.479 GeV , = 0.199. (12) 



Note that is positive, contrary to the other nonets 20]. The numerical values for the other constants k^, js 
and 75 g (which should be small if the large N c limit and chiral symmetry still applies approximately for the scalar 
nonet) are determined by a fit to data. 

Finally we discuss the mass relation in the octet sector. From the first three equations in IjlOfl we have: 

3M| g = 4M|, - M 2 ao + 3 7 s 8 Ml , (13) 

where the term proportional to 7s 8 includes deviations form the GMO limit. 

Glueball-quarkonia mixing is usually set up in the basis of fields N = y / l/2(uu + dd) = nn and S = ss instead of 
the octet and singlet fields So and Sg. The connection is given by 

So = y/2/3N + y/T/SS, S s = y/l/3N - y/2/3S. (14) 

Therefore, the scalar-isoscalar part of the Lagrangian © involving the fields N, S and G is rewritten as: 

C s = + M|)$ (15) 
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with 

M 2 N y/2f e 

V2f M% f | . (16) 
s f M 2 S 

The bare masses Mn, Ms and the flavor mixing parameter e are determined by inserting (|14f) into © with: 

M 2 N = Im 2 So + Im 2 S8 + ^zs. 



3 oo 3 o s 3 



Y(Ml-Ml s )-^ s . (17) 



The inverted relations are: 



\/2 , , 2 N 1 

^5 = IM* N +-Ml+--^:. 



- ! 2 1 , 2^2 
3 N 3 s 3 

So 3 JV I" 3 5 3 i 



igian CJ. 

BBBB0. 



*s = \/f (Mj-M|)-ie. (18) 

As we will show in the next section, as a result of the fit we will determine the values Mn, Ms, and e, from which we 
can deduce Ms , Ms s and zs or, equivalently by (|llf> . the parameters fc^ , 7s and 7s 8 of the Lagrangiaii 

The parameter / is the quarkonia-glueball mixing strength, analogous to the parameter z of Refs. 
The mixing strength z refers to the quantum mechanical case, where the mass matrix is linear in the bare mass terms. 
The connection between / and z, discussed in Refs. , leads to the approximate relation / ~ IzMg ■ 

Theparameter e induces a direct flavor mixing between the quarkonia states N and S. This effect is neglected 
in |a ISlul llfll l) where flavor mixing is considered to be of higher-order. However, a substantial N-S mixing in the 
scalar sector is the starting point of the analysis of Refs. [3 13- The origin of quarkonia flavor mixing is, according 
to connected to instantons as in the pseudoscalar channel, but with opposite sign (see also |3)- Such a phase 

structure is also found in the NJL model including the six-point t' Hooft interaction term full3)- The mixed physical 
fields are predicted to be a higher lying state of flavor structure — S]/y/3 and a lower one with [N + 5* a/2] / V3. 

Here we study the case e 7^ 0, more precisely e > 0, which leads to the same phase structure as in Ref. [3 13; but 
the quantitative results and interpretation will differ. 

For the glueball-quarkonia mixing we work in the limit of flavor blindness. The issue of flavor blindness breaking 
at the mixing level has been considered in 

BSE 

. This effect can be taken into account when introducing in the 
mass mixing matrix Ml, defined in (|15f) . an additional parameter r as: 

(Ml V2f-re \ 
Ml = y/2f ■ r Mq f . (19) 
\e f Ml) 

For r = 1 we regain the original expression, r 7^ 1 takes into account a possible deviation from this limit. Determination 
of this parameter on the Lattice ?J results in r = 1.20 ± 0.07, in the fit of 8] a value of r = 1 ± 0.3 is obtained. In 
the microscopic quark/gluon model of the value r ~ 1.1-1.2 is deduced. All these findings point to possible small 
deviations from the flavor-blind mixing configuration. Hence, in the following we will restrict to the limit r — 1. 
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The orthogonal physical states assigned as /i = /o(1370) , /a = /o(1500) and / 3 = / (1710) resulting from Eq. i(T5|) 
are obtained by diagonalization of M\ (Eq. l|16[ll with the transformation matrix B as 

BM 2 B T = M) = | " M| ] , (20) 

where the eigenvalues of M 2 represent the masses of the physical states /i = /o(1370) , /a = /o(1500) and 
/ 3 = / (1710) . The physical states \i), with i = fi, /2, /3, are then given in terms of the bare states as 

H)= B v\j)- ( 21 ) 

j=N,G,S 

In a covariant framework, where mixing of bound states is studied on a elementary level, it is not possible to define an 
orthogonal mixing matrix B [B ■ B T = B T ■ B = I3) connecting the physical to the bare fields [la [31L . However, 
as shown in the model of for the glueball-quarkonia system, deviations from orthogonality of the mixing matrix 
B are small, therefore justifying a Klein-Gordon mixing scenario as in the present approach. 

C. Strong and electromagnetic decays of scalar states 

The generic expression for the strong decay width of a scalar state s (both quarkonia and gluonium) into two 
pseudoscalar mesons p\ and P2 is given by: 

A 1 /2(M S 2 ,M„ 2 ,M„ 2 ) 

where X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2xz is the Kallen triangle function. The factor 7 spiP2 = 1/2 or 1 stands 
for identical or different particles in the final state. In the case of a degenerate isomultiplet an average over the isospin 
configurations is understood. 

The matrix elements M s ^ PlP2 are expressed in terms of parameters c^, c s m , c 9 d and c 9 m (see Eq. J2J0- The parameters 
c s d and c s m refer to the scalar quarkonia decays, c 9 d and to the direct glueball decays. The complete expressions for 
the two-pseudoscalar decay widths (matrix elements included) of the scalar resonances are given in the Appendix IbI 

The decay of the bare glueball states embedded in the physical /o states can proceed in two ways (see Introduction) . 
Mixing expressed by the parameter / corresponds to the conversion of a bare glueball to a quarkonia state, which in 
turn decays into a pseudoscalar meson pair. For / 7^ no physical /o state is a pure glueball, and the decays of the 
quarkonia components are driven by the amplitudes BiN and Bi$ 1 which depend on /. 

The direct decay of the glueball component, without proceeding via an intermediate scalar quarkonium state, is 
contained in the parameters c 9 d and c 9 m . For a microscopic description of this mechanisms we refer to Refs. \(i 
where the possible transition of the glueball to two pseudoscalar mesons is described by processes containing four 
internal quark/ antiquark lines. The parameters c d and c 9 m are attached to the gluonic amount BiQ of the state i, 
where i — /1, /a, /s< With the normalization adopted in Eq. (J2J) the limit c d — c 9 d , c s m — c 9 m refers to a direct glueball 
decay strength equivalent to the decay of a flavor-singlet quarkonia state. In the large N c limit the quarkonia decay 
constants c s d and c s m scale as N^ 2 , while the glueball decay constants c 9 d and c 9 m scale as N c [3o|. Large iV c arguments 
therefore suggest that the direct glueball decay is suppressed with respect to quarkonia decays. When performing in 
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the following a fit to data, we first study the large N c limit by setting the direct glueball decay parameters c 9 d and c 9 m 
to zero. 

The matrix element for the two-photon decay of the scalar (0 ++ ) state has the manifestly gauge-invariant form 

M s ^ 77 = e 2 g sll (g^ q x q 2 - q\q%) fyfai) e v (q 2 ) (23) 

where qi and q 2 are the photon four-momenta, p = qi + q 2 is the scalar state momentum and g sll is the S77 coupling 
constant. The decay width of the transition s — > 77 is given by 

Ts ^ = 32TmT £ \ M ^f\ 2 = \ a " 9%, M! , (24) 

pol 

where a = e 2 /(4ir) = 1/137 is the fine structure constant. 

The coupling constants <? S77 for the bare states So, Sg, G and for the isovector state Oq are directly calculated using 
the effective Lagrangian ©: 

9s n = ^/^ Ce ' 9s ^ = 3V3 Ce ' 
_32_ g __ 8 /- s 

5G77 — 3y / 3 Ce ' ^ a ° 7 ' r _ 6 

The parameter cj! refers to the quarkonia components, while cf contains the direct coupling of the glueball component 
to the electromagnetic fields. Latter coupling constant cf is supposed to be suppressed, since gluons do not couple 
directly to the photon field. However, an intermediate state of two vector mesons for example can in the framework 
of vector meson dominance generate a coupling of the glueball to the two-photon final state; this coupling is supposed 
to be suppressed and will not be considered in the numerical analysis. Using the identities for the field transforma- 
tions l|14|l and (|21|l we can derive the couplings for the bare N = yjl/2(uu + dd) and S = ss states and finally for 
the three physical scalar-isoscalar states i = fi, f 2 , and f& 

5 40 

V C e , 3i 77 — 

j=N,G,S 



9njj = 5s 77 - -g- V2 c s e , 3i 77 = ^2 Bij gjyj . (26) 



III. PHENOMENOLOGICAL FIT WITHOUT DIRECT GLUEBALL DECAY 

A. General considerations 

In the following we determine a best fit of the parameters entering in Eqs. J2J) to the experimental averages of 
masses and decay modes listed in Ref. [J. We first analyze the case of a non-decaying glueball, i.e. c? d — c 9 m = 0, 
where the decays are dominated by the quarkonia components (as in the original work of 6]) in line with large Nc 
arguments. The phenomenological analysis of Ref. ^(J confirmed this trend, but, as already mentioned, the fit of Q 
shows a strong contribution from the direct decays of the glueball configuration. 

The parameters of the model entering in the fit are the three bare masses -Mat, Mg, Ms, the two mixing parameters 
/ and e and the two quarkonia decay parameters c s d and c s m : 

M N , Mq, M s , f, e, c% c s m . (27) 

As an experimental input we use the following accepted values from Q: 
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(a) The scalar-isoscalar /o masses with the corresponding values of 



M /1S/0(13 70) = 1.35 ±0.15 GeV, (28) 
M> a =/ (i500) = 1.507 ± 0.005 GeV, (29) 
Mf 3 = M n W ) = 1.714 ±0.005 GeV. (30) 



(b) The partial decay widths of fa = /o(1500) with: 



T h ^ = 0.0380 ±0.0050 GeV, (31) 
r h ^KK = 0.0094 ±0.0017 GeV, (32) 
r/ 2 _ w = 0.0056 ±0.0014 GeV. (33) 

(c) The two accepted ratios for f% = /o(1710): 

Th-^fTh^EK = 0.20±0.06, (34) 
T f3 ^ vv /T h ^ K = 0.48 ±0.15. (35) 

r-i 

(d) The state /o(1710) has only been observed in the decays into two pseudoscalar mesons |l|. The decay into the 
final state 4ir, which can be fed by higher meson resonances, is suppressed :33|. We therefore impose the additional 
condition that the sum of partial decay widths into two pseudoscalar mesons (Tf 3 ) 2p saturates the total width (r/ 3 ) tot : 

(r /3 ) 2p = (r /3 ) fot = 140 ± 10 MeV. (36) 

Such a constraint is necessary to obtain meaningful total decay widths: without this condition on the full width a 
minimum for x 2 is obtained where (Tf 3 ) to t is larger than 1 GeV, a clearly unacceptable solution. 

(e) The two decay ratios of the 1=1 state ao(1450): 

Ta °^' = 0.35 ±0.16, (37) 

-L ciQ—f-Trri 



0.88 ± 0.23. (38) 



The total width of ao(1450) into two pseudoscalars is not known, because of the uncertainty of other decay widths as 

for LdTTTT . 

(f) The 1 = 1/2 state -K"q(145Q) decays dominantly into Kir with the corresponding width 

Tks->k* = 273 ± 51 MeV . (39) 

The only accepted average not included in the fit is Y f 2 ^ m * for the reason that the decay channel rjrj' is produced 
at threshold. Therefore, a significant distortion due to the finite width of the state is expected. 
For the N = 12 experimental values listed in Eqs. (J2SJ- l|39|) we perform a x 2 fit with 

2 



N=12 



X 2 = x 2 [M N , M G , M s , /, e, c s d , c s J = ^ 



i=i 



j^theory ^exp 
AA; 



(40) 
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where A° xp represents the i-th experimental result, AAi is its error and Al heo7y is the corresponding theoretical 
expression, depending on the parameters of 1)27(1 . 

For the case, where the direct glueball is suppressed (i.e. we set c 9 d = c 9 m = 0), two local minima for x 2 are obtained, 
the consequences of which we will describe in the following. The first solution was already analyzed in jl3| , while the 
second one is a novel solution with some peculiar characteristics. 

B. First solution and implications 



Fit results: From the first solution (I) we extract following fit parameters: 

M N = 1.455 GeV, M G = 1.490 GeV, M s = 1.697 GeV, 
/ = 0.065 GeV 2 , e = 0.211 GeV 2 , 
c d = 8.48 MeV, c s m = 2.59 MeV; xL = 29 - 01 ■ 
The corresponding fit results are reported in Table 1. 

Table 1. Fitted mass and decay properties of scalar mesons. 



(41) 



Quantity 


Exp 


Theory 


xl 


M h (MeV) 


1350 ±150 


1417 


0.202 


M h (MeV) 


1507 ±5 


1507 


~ 


M f3 (MeV) 


1714 ±5 


1714 


0.003 


T h ^ (MeV) 


38.0 ±4.6 


38.52 


0.011 


r /a ^s* (MeV) 


9.4 ± 1.7 


10.36 


0.322 


T h ^ m (MeV) 


5.6 ± 1.3 


1.90 


8.109 


Ffa-nriT/Tf 3 -iKK 


0.20 ±0.06 


0.212 


0.036 


r/ 3 -> w /Tj 3 _>kk 


0.48 ±0.15 


0.249 


2.446 




0.88 ±0.23 


0.838 


0.032 


TflQ — >n7]' /Tao — >7T7] 


0.35 ±0.16 


0.288 


0.150 


r K ,^ K „ (MeV) 


273 ± 51 


59.10 


17.590 


(r /3 ) 2P (MeV) 


140 ± 10 


143.27 


0.110 


Xtot 






29.01 



Bare masses: The bare non-strange quarkonia field N has a mass of Mn — 1.455 GeV, which is, as desired, similar 
to the scale set by the isotriplet combination ao(1450) with a mass of Mn^= 1.474 ± 0.019 GeV [l|. The mass of the 
bare glueball Mq = 1.490 GeV is in agreement with the lattice results [2| and with the phenomenological analyses 
of y,Q,^]]. The bare state S has a mass of A/5 = 1.697 GeV, which is about ~ 200 MeV heavier than the N state, 
an acceptable mass difference like in the tensor meson nonet. 

Mixing parameters For the glueball-quarkonia mixing parameter we get / = 0.065 GeV 2 , which by the approximate 
relation / ~ 2zM G 0, 12 Lcorresponds to z ~ 21.8 MeV. 

The results of Refs. J 3 E3 are z = 85 ± 10 MeV, z = 80 MeV and z ~ 62 MeV, respectively, i.e. of the same 
order, but larger. The introduction of additional flavor mixing between the quarkonia configurations in the fit, as 
done here, leads to a reduction of the strength parameter /. The lattice result of Ref. 7] with 43 ± 31 MeV is in 
agreement with the present evaluation, but has a large uncertainty. A mixing strength of the same order is found in 
the lattice evaluation of [l^ . 
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The flavor mixing parameter resulting from the fit is e = 0.211 GeV 2 . In the limit / = the mixed physical states 
are = 0.97 \N) + 0.26 \S) and |/ 3 ) = -0.26 L/VV+ 0.97 \S) (and, of course, |/ 2 ) = \G)). The phase structure of the 
mixed states is, as discussed previously, as in |23l I28M 29). But here the strength of flavor mixing is smaller, resulting 
in mixed states, which are dominantly N or S. The influence however of (an even small) flavor mixing in strong and 
electromagnetic decays may be non- negligible. 

Mixing matrix: The mixing matrix B relating the physical to the bare states in the present fit is expressed as: 

/ 1/0 = |/o(1370)) \ / 0.86 0.45 0.24 \ f \N) = \nn) \ 

|/a) = |/ (1500)> = -0.45 0.89 -0.06 \G) = \gg) . (42) 
V l/a) = |/ (1710)> / V -°- 24 -°- 06 °- 97 / V I 5 ) = \ Ss ) J 

The physical resonances are dominated by the diagonal bare components, qualitatively in line with Refs. 
Since the glueball does not contribute to the decay, the relative phase with respect to the quarkonia components is 
at this stage irrelevant. By inverting / — > — / we would find the same results for the decays, but opposite glueball- 
quarkonia phases. In turn, the relative phases of the N and S components are not symmetric under e — ► — e. As 
discussed above, in /o(1370) they are in phase, while in /o(1710) they are out of phase. The state |/o(1500)) behaves 
like a N state with a decreased width, while the S component is small ^1 13 . Thus the decay into KK is smaller 
than for the tttt channel. In the present solution (I) the N and S state components are in phase contrary to the results 
of . However, the other solution (II), presented later on, shows again an opposite phase in / 2 = /o(1500), but 

with a large ss amount. 

Large N c constants: From the fit parameters of (|41|) we determine 75,,, 7s 8 and by using l)18|l and I jllf l: 



1.479 GeV, = 0.199, 
0.225, 7s - 0.236, k s m 



-0.818 



(43) 



The values of 75,, and are smaller than in 20]. Also, the violation of the GMO relation, encoded in 7^, is small, 
indicating that higher order corrections in the chiral expansion are possibly not too large to invalidate the present 
study. In the present scenario the smallness of the glueball mixing parameter / can also be interpreted as a small 
violation of the large iV c limit. The present results show that large N c and chiral symmetry can, although violated at 
some level, be a useful guideline to scalar meson physics. 

Resonance fn( 1370): The experimental uncertainties of the /o(1370) resonance are large, no average or fit is 
presented in [l|. The main problem connected with this resonance is its large width (200-500 MeV) and its partial 
overlap with the broad low- lying a = /o(400 — 1200). However, the results from WA102 [3^ indicate a large N = Tin 
component in its wave function. Results from Crystal Barrel (summarized in |2| and subsequently analyzed in [3^^ 
confirm such a trend (see also for a recent review). 

Table 2. Decays of fx = / (1370). 



Quantity 


Exp (WA102) 


Theory 


(r*)^ (MeV) 


0.46 ±0.19 
0.16 ±0.07 
"small" 


0.34 
0.06 
166 
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Predictions for the two-pseudoscalar decay modes are in acceptable agreement with the results of WA102 as shown 
in Table 2. The measured ratio r/ 1 _ > 4 7r /r/ 1 _ >lr7r = 34.0^g 2 [3], although the errors are large, points to a dominant 
4-7T contribution to the total width. Our prediction gives however a sizable contribution of the two-pseudoscalar decay 
mode and is therefore not in agreement with such a large 47r decay mode. 

In the original work of Q, a quarkonium Tin state has a very large two-pseudoscalar width with T^n-nnr = 270 ± 
25 MeV, T- n _^ K = 195 ± 20 MeV and T- m ^ m = 95 ± 10 MeV, i.e. {T- m ) 2P ~ 500 MeV. The /i = / (1370) 
is in 6] dominantly nn, therefore one expects a large value for (Tf 1 ) 2P , in contrast to the presented experimental 
analyses listed above. A large value for (Ynn)2P is also predicted in [10(: for the mixed state /o(1370) one has 
(Tf 1 ) 2p — 115.7 MeV comparable to the present study. 

On the contrary in the study of 8], small two-pseudoscalar partial widths are obtained by the following mech- 
anism: the glueball decay amplitude in /o(1370) to two pseudoscalar is large and interferes destructively with the 
Tin component (the phases quarkonia-glueball in |8j are inverted with respect to l|42|l). As a result |8j the two- 
pseudoscalar decay width (Tf 1 )2p is smaller than (Tf 2 )2p with (T f 2 )2 ps /(T f 1 )2ps = 10.0 ± 3.0. At the same time 
(^hhp/(T fl )2P = 0.7 ±0.2 is obtained. 

In the present fit with an inert glueball (as in the original work of 6] and as in [10| , where the glueball is allowed to 
decay, but the quarkonia components still dominate) we find the following decay widths into two pseudoscalar pairs: 

(T fl ) 2P = 166 MeV > (r /3 ) 2P = 143MeV > (r /2 ) 2P = 51 MeV. (44) 

An analysis by Crystal Barrel also indicates sizable partial decay widths of the Att decay channels: _,.<,■<, = 
120.5 ± 65 MeV and Tf 1 ^ rpp = 62.2 ± 28.8 MeV. The same analysis gives the following two-pseudoscalar partial 



widths 36, 3 



ft— >7T7T 



21.7±9.9MeV, T 



fi—>KK 



(7.9±2.7MeV) to (21.2±7.2 MeV), r /x _ 



0.41±0.27MeV. 



On the contrary, the analysis of |38i] reports the ratio r/ 1 ^ 7r7r /(r/ 1 ) tot = 0.26±0.09, pointing to a large hit (ergo to a 
large two-pseudoscalar) partial decay width for /o(1370). This experimental result is therefore in disagreement with 
the analysis of |36j. New results on /o(1370) would be crucial to disentangle the scalar puzzle and to understand if a 
destructive glueball/quarkonia interference as in v8| is necessary. 

Resonance /o(1500): The theoretical partial widths of the /o(1500) are in good agreement with the data (see 
Table 1) apart from a slight underestimate of the 2n channel. We also obtain T f 2 ^ vri > ~ 0.036 MeV as compared to 
the experimental value of r ' f 2 ^m' = 2 ± 1 MeV. Taking into account the finite width of the resonance will lead to an 
increase of the theoretical value. 

Resonance /o(1710): For the decays of /o(1710) we summarize our results compared to the data of WA102 33] in 
Table 3. 



Table 3. Decays of / 3 = / (1710). 



Quantity 


Exp (WA102) 


Theory 


r , / 3 - > ^Rr/r/ 3 -„r7r 


5.0 ±0.7 


4.70 


F/s— >rw/Ff 3 — >7nr 


2.4 ±0.6 


1.17 


F/a ->rrn' /r/ 3 -^-ktt 


< 0.18 


1.59 


(T h ) 2P (MeV) 


" dominant" 


143.27 



The first two ratios, already included in the fit of Table 2, can be reproduced. The theoretical ratio T/ 3 ^^^' /Tj 3 ^^ w , 
which is not included in is in complete disagreement with the WA102 result. The dominance of the r\r{ mode over 
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7T7T is a solid prediction in the framework of solution I, which does not depend very much on the choice of parameters. 
A confirmation of the experimental result could possibly hint at a different mixing scenario or at a sizable role of 
direct glueball decay. 

Resonance ao(1450): The ratios of two-pseudoscalar decay modes of ao(1450), included in the fit of Table 2, are 
well reproduced. The prediction for the two-pseudoscalar width of (T ao )2P = 84.26 MeV is smaller than the total 
width of 265 ± 13 MeV. However, the experimental ratio (r ao ^ W7r7r /r ao ^ 7r ^) is not known: no average or fit is listed 
in 111. The experimental value from [3JJ, which is 10.7 ± 2.3, would imply a dominant ojtttt mode and in turn a 

, nn n 

rather small two-pseudoscalar partial decay widths. This finding is in disagreement with the results of [6J,|40J. In |6J 
a value of (r ao ) 2 p = 390 ± 110 MeV is found, in the work of Q one has (T ao ) 2P = 420 - 940 MeV for a 
masses in the range of 1200 — 1400 MeV. For our value for (T ao )2p we obtain the estimate: T l o->w7r7r /r a — >7r?j = 

[(Fao)tot ~ (r oo )2P]/ Fao-K^ ~ 4.5 . 

Resonance -?lq(1430): Our result for Tk^k-k underestimates the experimental value by a factor of about 5 (see 
Table 1). Furthermore, for the additional Kij decay channel we get T k*^,Kt]/^ K'^irK — 0.026. 

In |4pJ a value of Tk*->k-k = 340 MeV is predicted, but, as discussed above, (r ao ) 2P is of the order of 1 GeV, much 
larger than the full width. Similarly, in [6J with T K *^ Ki: = 200 ± 20 and (r Qo ) 2 p = 390 ± 110 MeV the first result 
underestimates while the second overshoots the experimental value. A full analysis in the 3 Pq model |4l| results in 
Fk*^ttK — 166 MeV and Tn^^^ — 271 MeV; unfortunately the resonance ao(1450) is not discussed in 41]. The 
authors of |4l| also tried to adjust Tk^ttK to its experimental value, and then calculate the 2n partial width of a 
N state, obtaining IV— ~ 450 MeV. The last result implies a very large two-pseudoscalar and full width for a N 
state. A full experimental determination of all relevant decay modes involving ao(1450) and if □ (1430) would certainly 
help to clarify this issue. We refer to section |TlI Dl for a further discussion of this problem. 

If /o(1370) is dominantly nn, as in P. R. IH. Ifll llfl ll^ . there is, as discussed above, an incompatibility of the present 
experimental small two-pseudoscalar partial decay widths and various model calculations. At the same time, a 
consistent understanding of the isodoublet states i4Tg(1450) and the isovectors ao(1450) is still incomplete. 

Two-photon decays: As a further consequence we discuss the two-photon decay rates of the scalar resonances. We 
assume that the coupling c 9 e is suppressed with respect to c s e , i.e. we set the glueball-photon coupling c| to zero. The 
ratios of radiative decay widths as a prediction of the fit are: 

r/^2 7 : r /2 ^ 27 : r /3 ^ 27 : r a <^ 27 = 1 : 0.305 : 0.002 : 0.471 , (45) 

which are independent of the coupling c*. The result for r/ 2 _ > 27/r/ 1 _ > 27 is in qualitative agreement with the results 
of . The ratio Tfo-^/Tfa-^, however, is considerably smaller than in the previous works. The suppression 

of rj3_»27 originates from the destructive interference between the N and S components, which in turn is traced 
to the flavor mixing with e > in accord with the phases of |3 |2 • Another interesting prediction is the ratio 
r a o_^2 7 /r/ 1 _ + 27, which is relatively large. 

The experimental status of the two-photon decays is still incomplete. For the /o(1370) two values are indicated 
in PDG2000 ^^|: 3.8 ± 1.5 keV and 5.4 ± 2.3 keV. However, it is not clear if the two-photon signal comes from the 
/o(1370) or from the high mass end of the broad /o(400 — 1200). The PDG currently dj seems to favor this last 
possibility, but the data could also be valid for the /o(1370). We therefore interpret the two experimental values as an 
upper limit for the two-photon decay width of the /o(1370). Signals for two-photon decays of /o(1500) and /o(1710) 



have not yet been seen, the following upper limits are reported 
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r/ (i5oo)^2 7 (r/ (i5oo)^Wr/ (i5oo)tot) < 0.46 keV, 
r/ (i7io)^2 7 (r /o (i7io)^KK/r /o (i7io) to t) < O.ll keV. (46) 

Using the known branching ratio T ^^qq^^tt f (i500)tot one gets 500)^27 < 1-4 keV 34]. An accepted fit for 
r / (i7io)^Ki?/ r /o(i7io)tot is not reported in [l|. Using the value from |4j| with r /o(mo) ^ Ki? /r /o(mo)tot = 0.38±o;°3 
we find an upper limit of the order of ^ f (mo)-f2j ~ 0.3 keV. 

For an absolute prediction of the two-photon decay widths we use c s e = 0.0138 GeV -1 as determined in the model 
approach of Ref. |l2( For the non-strange quarkonium state we get Tn^j — 0.969 keV, while for the isovector and 
mixed scalars we have: 

T fl ^ 2l = 0.703 keV, r /a ^ 27 = 0.235 keV , 

r /3 ^2 7 = 0.002 keV, r a o^ 27 = 0.362 keV . (47) 

The results for the mixed states are below the current upper limits (the original results presented in jl3] contain a 
slight missprint, the correct numbers are reported here. The physical considerations are not affected from this slight 
change) . 

The estimate for the 27 decay of the bare quarkonium state N = fin of 0.969 keV is smaller than the one of j^ . 
where the following expression has been used: 

r M7 (0 ++ ) = k ( ^Sy ) IW 27 (2++) (48) 

The coefficient k is 15/4 in a non-relativistic calculation, but becomes smaller when considering relativistic cor- 
rections 0| . I n H| a range of values for k from 2 to 15/4 is considered. Our chiral Lagrangian approach 
combined with [l2| points to a smaller value of k. Using our result for Fjv^2 7 (0 ++ ) and taking the value 
rAr^2 7 (2 ++ ) = 2.60 ± 0.24 keV jl| at M N {2++) = 1.27 GeV we get k ~ 0.25. This result is model dependent, 
since it relies on the parameters for the covariant description of the scalar mesons used in 12]. A fully covariant 
treatment may imply strong deviations from the non-relativistic limit. 

Discussion: The largest contribution to x 2 is due to the underestimate of the Kq — > Kir width: from Table 1 we 
have x% ot /N = 2.42. When excluding the data point for Tk*—>Kit in the fit, a very similar minimum compared to 141|) 
is found with: 

M N = 1.442 GeV , M G = 1.485 GeV , M s = 1.695 GeV , / = 0.080 GeV 2 , 

e = 0.225 GeV 2 , c s d = 8.12 MeV , c s m = 3.57 MeV ; Xtot = n - 19 ■ ( 49 ) 

In this case we have Xtot/N = 1.02, corresponding to a good description of data. The discussion about the isoscalar 
states and ao(1450) remains unchanged. 

The underestimate of Kq — > Kir constitutes an open problem of the scalar analysis. As already discussed, a 
consistent understanding of the complete scalar nonet is lacking in other approaches as well. We will further discuss 
this issue in section IlII Dl 

Aside from this difficulty, the rest of the accepted data in [j] is well described. The quality of the current fit 
seemingly excludes a sizable direct decay of the scalar glueball component. Concerning results for data, which are not 
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reported as average or fit in [lj], the situation is less clear: the predicted full two-pseudoscalar width of f\ = /o(1370) 
is larg e, when confronted with the WA102 result (see Table 2 and Refs. but also the different result 

of |38j). The ratio r , / 3 _>, 7 ,j//r'/ 3 _> 7r7r (Table 3) is also problematic in the present solution. An accepted average, in 
particular for these values, would help in clarifying these issues. 



C. Second solution and implications 



Fit results: Solution II is obtained for following fit parameters: 

M N = 1.298 GeV , M G = 1.513 GeV , M s = 1.593 GeV , / = 0.400 GeV 2 , 
£ = 0.015 GeV 2 , c s d = 7.48 MeV , c s m = 6.42 MeV ; X tot = 24.61 . 
The corresponding results are reported in Table 4. 

Table 4. Fitted mass and decay properties of scalar mesons. 



(50) 



Quantity 


Exp 


Theory 




M h (MeV) 


1350 ±150 


1142 


1.924 


M h (MeV) 


1507 ±5 


1508 


0.023 


M h (MeV) 


1714 ±5 


1713 


0.0254 


r A _™r (MeV) 


38.0 ±4.6 


37.31 


0.019 


Tf^HK ( M eV) 


9.4 ± 1.7 


10.08 


0.167 


T h ^ m (MeV) 


5.6 ± 1.3 


4.70 


0.477 




0.20 ±0.06 


0.216 


0.071 


^ *fa->vvfi f 3 —>KK 


0.48 ±0.15 


0.248 


2.400 




0.88 ±0.23 


1.078 


0.741 


Tao — 'TTf]' /"^ a o — 


0.35 ±0.16 


0.291 


0.134 


Tk^k, (MeV) 


273 ± 51 


53.75 


18.480 


(r /3 ) 2p (MeV) 


140 ±10 


143.94 


0.155 


Xtot 






24.61 



A first look to (|f>C)|) shows some peculiar differences when compared to the set of 141(1 . In the following we discuss the 
implications and the differences of this second solution. 

Bare masses: The bare masses Mjv and M$ are smaller than in solution I, their mass difference is still around 
200 GeV (Eq. (|50J) i. The bare glueball mass is about ~ 1.5 GeV, as before, but now is much closer to Ms- This small 
mass difference leads to a strong mixing between the glueball and the bare S =ss state. 

Mixing parameters: The quarkonia flavor mixing e = 0.015 GeV 2 is very small in this solution and has practically 
no influence on the phenomenology. On the contrary / = 0.400 GeV 2 is much larger, leading to a strong glueball- 
quarkonia mixing. In this respect there is a clear difference between the two solutions. Using the approximate 
relation / ~ 2zMq we find in this case z ~ 130 MeV, which is larger than the results from other works listed 

in Sec. III.B.3, but still in qualitative agreement. 

Mixing matrix: The mixing matrix B, relating the physical to the bare states, for the second solution is expressed 
as: 

/ I/O = |/ (1370)> \ 

|/ a ) = |/ (1500)> 
\ |/ 3 ) = |/ (1710)> J 



0.81 0.54 0.19 
-0.49 0.49 0.72 
-0.30 -0.68 0.67 



\N) = \nn) 

\G) = \gg) 
\S) = \ss) 



(51) 
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The large mixing parameter / causes the glueball configuration to be spread out among the /o states: /i = /o(1370) 
is still dominantly nn, fa = /o(1500) is mostly ss, but with a sizable out-of-phase nn amplitude (the opposite phase 
of nn and ss was first considered mechanism to explain the large irir/KK ratio). In = /o(1710) both 

the gluonium and the ss components are large, where the gluonic component is slightly larger. Remarkably, although 
the bare level ordering is still Mn < Mo < Ms, the largest ss amount is contained in fo, while the largest gluonic 
component is present in fa. The mixing matrix resembles some features of the results of |7j, although the bare level 
ordering is different. 

Large N c constants: From the parameters in l|5U|) we determine 7s , 7s 8 and by using i|ll|) and 11811 : 



M s = 1.479 GeV, = 0.199, 
7Ss = 0.139, 75o = -0.032, k s m 



-0.786. 



(52) 



The same considerations as for solution I also hold here. 

Resonance /o(1370): The results for /i = / (1370) are summarized in Table 5. 

Table 5. Decays of /i = / (1370). 



Quantity 


Exp (WA102) 


Theory 


(T fl ) 2p (MeV) 


0.46 ±0.19 
0.16 ±0.07 
"small" 


0.27 
0.02 
56.79 



The main difference with respect to solution I concerning /i = /o(1370) is the decreased theoretical two-pseudoscalar 
width (mostly caused by the smaller physical mass, see Table 4). In the present fit we find the following decay widths 
into two pseudoscalar pairs: 



(T h ) 2P = 144 MeV > (T fl ) 2P = 57 MeV > (T h ) 2P = 52 MeV. 



(53) 

Here the decay pattern is in better agreement with the analysis of |36j, |37j than the one of solution I. 

Resonance /o(1500): The theoretical partial widths of /o(1500) are in rather good agreement with the data (see 



Table 2). We also obtain r / 2 -* w ' = 1.5 MeV as compared to the experimental value of Tf 2 ^ nn > = 2± 1 MeV without 
invoking further threshold corrections. 

Resonance /o(1710): For the decays of /o(1710) we summarize our results compared to the data of WA102 
Table 6. 



Table 6. Decays of / 3 - / (1710). 



Quantity 


Exp (WA102) 


Theory 


f3 -+KK 1 r f3 -» 7T TT 


5.0 ±0.7 


4.63 


^f3— /r/ 3 ->7r7r 


2.4 ±0.6 


1.15 


^/a -*V*l' AVa -^ktt 


< 0.18 


0.36 


(r /3 ) 2p (MeV) 


" dominant" 


143.94 



The theoretical r/ri' /nn ratio is now smaller than in solution I. Although the prediction is still larger than the upper 
limit set by WA102, it does not represent such an evident mismatch. 
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Resonance ao(1450): The two experimental ratios are satisfactorily described. Otherwise, the discussion of solution I 
is still valid here. 

Resonance A"q(1430): The underestimate of Tk*^>k-k is also present in the second solution (see Table 4). Further- 
more, for the additional Krj decay channel we get T K*—>Kri/^ K*—>TtK — 0.05065. 
Two-photon decays: The two-photon decay ratios resulting from solution II read: 

r /l ^ 27 : r /2 ^ 27 : r /3 ^ 27 : r a o^ 27 = 1 : 0.253 : 0.055 : 0.493. (54) 

The two 7 results are in qualitative agreement with 8, 121. 



Following the arguments of the previous section by using the result of [l^ we find: 

= 0.350 keV, T h ^ 2y = 0.088 keV, 
r /3 ^2 7 = 0.019 keV, r Q o^2 7 = 0.172 keV. (55) 

As before, the rates are below the upper limits set by current data, but they depend on the model-dependent evaluation 

of Q. 

Discussion: If in this case we exclude the Kq — > Kit mode from the fit (which generates the by far largest 
contribution to x 2 ), the following set of parameters is found: 

M N = 1.313 GeV, M G = 1.511 GeV, M s = 1.594 GeV, / = 0.395 GeV 2 , 

e = 0.002 GeV 2 , c s d = 7.30 MeV, c s m = 6,61 MeV; xLt = 5 - 97 - ( 56 ) 

We then obtain Xtot/N = 0.54, corresponding also in this case to a good description of the remaining data. Again, 
there is no phcnomenological evidence for a direct decay of the glueball component. 

Apart from the decay Kq — ► Kir, both solutions (I and II) describe the data l|28|) - (|38[) well. There are, however, 
differences when comparing to the WA102 results. The second solution comes in this respect closer to these data, but 
the experimental results are not yet conclusive. 



D. _R- *(1430) 

The most striking mismatch with the data is in both analyzed scenarios I and II the theoretical underestimate of 
the Kq — > Kit mode (Tables 1 and 4). The corresponding partial Xk*—>k-k * s ^ ar ^ ne dominant contribution to 
Xtot- This mismatch also holds when including direct glueball decays in the analysis as shown in the next section. 

As a further attempt, following [4 1| . we can also pursue another strategy: in a first fit we fix the quarkonium decay 
parameters c s d and c s m in order to reproduce Kq — > Kir (Eq. (|39|l ) and the two ratios for ao given in (|37|l and l|38|l . 
Then we obtain c s d — 17.94 MeV (larger by a factor 2 when compared to the previous fit in Tables 1 and 4). and 
c s m = 7.35 MeV. With these values one has Y K *^K = 281 MeV, ^ aa ^K K /^ a a ^ v = 0.88 and T ao ^ vl /T ao ^ n = 0.30 
and r ao ^ 2 p = 381.17 MeV (already larger than the experimental result). 

The corresponding decay width of the N state into tw o p seudoscalars is Tjv^2P ~ 900 MeV (for a mass of 
Mjv ~ 1.4 GeV). We then find similar results as in [y, |4lJ, implying that the discussed trend is rather model 
independent. If at this stage we fix c s d = 17.94 MeV and c s m = 7.35 MeV and we make a second fit of the remaining 
free parameters (Mn, Mq, Ms, /, e, <A, c^J to the values reported in (|28|l - l|36|) . we find minima with Xtot ~ ^00, 



18 



which are clearly unacceptable. This is a further confirmation of the incompatibility of results in the scalar mesonic 
sector. If we attempt to reproduce the width of Kq(1430), the other results are off. As noted above, this problem 
arises in other models as well. 

In [J only the Kq — > Kir mode is quoted in the list of the decay modes. Our suggestion is that a strong coupling to 
the scalar mesons below 1 GeV, Kq (800) and a takes place. Both states are very broad and could eventually influence 
the decay strengths of ^(1430). The theoretical description of such a phenomenon is however beyond the goals of 
the present tree-level study. 

IV. PHENOMENOLOGY INVOLVING THE DIRECT GLUEBALL DECAY MECHANISM 

A. Flavor symmetry limit 

The analysis of the previous section did not reveal a phenomenological need to include the direct two-pseudoscalar 
decay of the glueball component. In this section we analyze this additional mechanism when including it in the fit. 
We start first by including the interaction term proportional to c d , while neglecting the flavor symmetry breaking 
contribution c 9 m . 

The ratio c 9 d /c s d is a measure of the direct glueball decay strength. The analogous quantity used in Q is, because of 
the different normalization, related as c 9 d /c s d < — ► y/3/2r2- For the different solutions in the fit of Q the parameter ri 
varies between 1 and 5. Note that the case Ti = 1 corresponds to a direct glueball decay strength into two pseudoscalars 
~ 1.22 larger than the quarkonium strength. In 13] the bare masses have been kept fixed when the glueball decay 
parameters have been introduced. Here we release this constraint by leaving the bare masses free. In the fit we again 
find two different solutions, which correspond to the ones (I and II) analyzed in the previous section. 

The x 2 minimum corresponding to solution I but now with inclusion of c d in the fit is: 

M N = 1.416 GeV, M G = 1.493 GeV, M s = 1.694 GeV, 
/ = 0.075 GeV 2 , e — 0.241 GeV 2 , c s d = 8.73 MeV , 
c s m = 1.48 MeV, c d = -0.94 MeV; X L - 28.70 . (57) 

A comparison with the previous fit l|41(l shows that the parameters practically did not change: with c d = —0.94 MeV 
the direct glueball decay is suppressed, resulting in |c?/c^| = 0.11 1 . Also, the total \ 2 is only slightly smaller than 
in (HJ. 

The minimum corresponding to solution II reads: 

M N = 1.401 GeV, M G = 1.461 GeV, M s = 1.609 GeV, 
/ = 0.376 GeV 2 , e = -0.082 GeV 2 , c s d = 7.63 MeV, 
c s m = 7.09 MeV, c 9 d = 1.82 MeV; xlt = 22.87 . (58) 

A comparison with H50[) shows a clear similarity, although the bare masses are somewhat shifted. Again, direct glueball 
decay is suppressed with |c^/c^| = 0.24. 

The inclusion of c d in the fit does not lead to a drastic change of the previous fit parameters. The bare glueball 
decay strength is strongly suppressed with respect to the quarkonium one, thus in agreement with the analysis of the 
previous section and with large N c considerations. 
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B. Flavor symmetry breaking and lattice results 

Inclusion of the flavor symmetry breaking term in direct glueball decay, that is the term with c^, results in nine free 
parameters (the set of ()27[) and c 9 d , c^J. A direct fit of these parameters to the data generates various minima, which 
are not well pronounced. For example, solutions are found where the flavor symmetry breaking term in the glueball 
sector c 9 m is exceedingly large, dominating the decay mechanism. Instead, to study the effect of flavor symmetry 
brea king in glueball decay we resort to a first lattice study 24] to fix the decay parameters c 9 d and c 9 n . 

In [24J a full two-pseudoscalar decay width of the scalar glueball with about 100 MeV is deduced. The corresponding 
mass of 1.7 GeV led the authors of [24| to interpret the resonance /o(1710) as mainly gluonic. However, in the cited 
lattice analysis it is not clear if the glueball decay mechanism occurs partially by mixing with scalar quarkonia 



(here parametrized by /), or by direct decay (parametrized by c 9 d and c^J. In jjj a scenario is studied, where the 
lattice results of j^l are explained by mixing only (hence c 9 d and c 9 m are set to zero). The physical state /o(1710) is 
mainly gluonic, but because of mixing, it acquires a large ss amount, which it turn explains the decay pattern. The 
corresponding mixing matrix is then similar to the results of Q (and to some aspects of our solution II). 

In it is stated that the amplitudes deduced in [2^ probably include significant contributions from mixing of the 
scalar glueball with quarkonium, although not proven. If this is the case, only the mixing mechanism contributes to 
glueball decay and then we are back to the previous solutions, where the constants c d and c 9 m are negligible. Here 
we also intend to investigate the opposite case, where the decay couplings calculated in [2^] arise from direct glueball 
decay. 

The decay widths of the glueball can be expressed as Q (see also Q): 

3AV 2 (M§,M^,M|) 2 

Vg — 32^m| iVo^Mp} , 

L G^KK - 8?rM 3 WGKK^PS > 

where M p = 775 MeV is the p mass. The lattice results for the decay constants yc-m*-, Vgkk anc ^ VGvv are: 

y G ^ = 0.834™, y Gl?K = 2.654l8:tg, V(*m = 3-099±g:^|. (60) 

Note that in the flavor-symmetry limit we would expect dg-k-k ■ Vgkk ■ VOvv = 1:1:1. Although the errors are 
large, the lattice results show a sizable deviation from this limit. As already noted, it is however not clear, if and to 
what extent mixing with quarkonia is included in these amplitudes. Interpreting the lattice results in the context of 
the direct glueball decay mechanism implies a large symmetry violation parameter c 9 m . 

The corresponding decays can be derived from the expressions of Appendix B by setting the glueball-quarkonia 
mixing to zero, that is by considering the scalar-isoscalar decay for Mj = Mq and BiN — Bis = and, of course, 
BiG = 1. The explicit expressions for the tttt and KK modes of the direct glueball decay are: 

Tg — - 32^Ml \-^([ M G-2M 7I \c d + 2M^ m )^ , 

X^{M G ,Ml,Ml) \ 2 _ V 

SnM G \V3F2 {[ G 2M K \c d + 2M K c m )j 



>-„• = r;;r ,JK ' <i-^(i^ - -^i <■:, f^d}' (ei) 



20 

Matching the expressions for the tttt and KK decay modes of (|61fl to (I59|l using Mq — 1.7 GeV (as in j^J]) we obtain 
for the decay constants: 

c 9 d = 1.34 MeV and c g m = 24.6 MeV. (62) 

For a bare glueball mass of Mq = 1-5 GeV we find rather similar values of c d ~ 1.72 MeV and c 9 m = 25 MeV; hence 
we have a rather slight dependence on Mq within a reasonable range of values. In the following we take the values 
evaluated at Mq = 1.7 GeV. 

Using the values of Ij62(l we can also determine the r\r\ decay amplitude. Compared to the lattice result of 
UGirq M p = 2.40lg'33 MeV we get the value of 2.025 GeV. The corresponding decay widths for the bare glue- 
ball are (for Mq = 1.7 GeV): 

ro-mr = 7.23 MeV, r G _^ x = 80.61 MeV, (63) 
r G ^„ = 18.35 MeV, T G ^ W = 11.73 MeV. (64) 

For a lower value of mass Mq the only significantly affected mode is 7777', since threshold effects become important 
(note that the 7777' mode is entirely generated by the flavor symmetry breaking term proportional to c^J. 

With the direct glueball decay including flavor symmetry violation fixed by l|62l) . we now perform a fit with the 
remaining free parameters Mjy, Mq, Ms, /, s, c s d , c s m . Two solutions (III and IV) arc obtained, which correspond to 
the bare level orderings Mjy < Mq < Ms and M/y < Mq < Ms, which we analyze in the following. Although the 
direct glueball decay dominates in these cases, the two solutions have similarities to the ones discussed in detail in 
the previous section. 

C. Third solution and implications 

Solution III is obtained for the set of parameters: 

M N = 1.359 GeV, M G = 1.435 GeV, M s = 1.686 GeV, 
/ = 0.212 GeV 2 , e — 0.277 GeV 2 , 

c d = 8.28 MeV, c s m = 7.21 MeV; xlt = 21-56 , (65) 
where the fit results are listed in Table 7. The mixing matrix is similar to the one of solution I and explicitly reads: 



= |/ (1370)> \ 
|/ a > = |/ (1500)) 
l/a) = |/o(1710)> / 



0.79 0.56 0.26 
= I -0.58 0.81 0.02 
-0.20 -0.16 0.97 



/ \N) ee |nn) 

\G) = \gg) 

V \S) EE \§ s ) 



(66) 



The decay rates for fi — /o(1370) and fa — /o(1710) as compared to the WA102 data are summarized in Tables 8 
and 9. 

The two-photon decay ratios resulting from solution III read: 



r /l ^ 27 : r /2 _ 27 : r /3 ^ 27 : r Q o_ 27 = 1 : 1.018 : 0.025 : 0.494. 



(67) 
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The results of the fit summarized in Table 7 are acceptable, apart from the already discussed underestimate of the 
K* decay width. The corresponding prediction for the WA102 data on /o(1370) has problems: the predicted ratio 
r f t _^K k /T /i— >7T7r i s larger than unity and the full two-pseudoscalar decay width is very large (we refer to the discussion 
of solution I on the issue of the latter point). For the state /o(1710) we obtain a large ratio T f 3 ^ vv > /T f^^^ , again 
as in solution I, in contrast to the WA102 result. For the two-photon decays we have a large ratio T f 2 ^2-y 

Table 7. Fitted mass and decay properties of scalar mesons. 



Quantity 


Exp 


1 heory 


2 

Xi 


71 /r /TV X T T\ 

Mf t (MeV) 


1350 ±150 


1242 


0.519 


M h (MeV) 


1507 ±5 


1507 


0.003 


M f3 (McV) 


1714 ±5 


1714 


0.011 


T f2 ^ (MeV) 


38.0 ±4.6 


38.50 


0.010 


r /a ^KK ( Mc v) 


9.4 ± 1.7 


10.38 


0.332 


Y h ^ m (MeV) 


5.6 ± 1.3 


3.65 


2.252 




0.20 ±0.06 


0.197 


0.002 


^fs—WI /Tf 3 -*KK 


0.48 ±0.15 


0.314 


1.221 




0.88 ±0.23 


1.079 


0.745 


Tao — /Fao — >7rr] 


0.35 ±0.16 


0.291 


0.134 


Tk^k, (MeV) 


273 ±51 


71.02 


16.221 


(r /3 ) 2p (MeV) 


140 ± 10 


143.3 


0.109 


Xtot 






21.560 



Table 8. Decays of fi = /o(1370). 



Quantity 


Exp (WA102) 


Theory 


(T fl ) 2p (MeV) 


0.46 ±0.19 
0.16 ±0.07 
"small" 


1.10 
0.17 
193.5 


Table 9. Decays of f 3 = / (1710). 


Quantity 


Exp (WA102) 


Theory 


^ fa -*K K 1 1" h -» kit 

^/3— *7]T] f 3 ^TVK 

r/ 3 — trm'/Tfa—HTn 

(r /3 ) 2P (MeV) 


5.0 ±0.7 
2.4 ±0.6 
< 0.18 
"dominant" 


5.08 
1.59 
2.01 
143.3 



D. Fourth solution and implications 

an inverted bare level ordering and a small glueball-quarkonia mixing with 

= 1.392 GeV, M G = 1.712 GeV, M s = 1.452 GeV, 
= -0.050 GeV 2 , e = 0.232 GeV 2 , 

= 6.66 McV, c s m = 5.84 McV; Xtot = 26.330 . (68) 



The solution IV corresponds to 

M N 
f 

4 
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Table 10. Fitted mass and decay properties of scalar mesons. 



Quantity 


Exp 


Theory 


Xi 


M h (MeV) 


1350 ±150 


1330 


0.017 


Mf 2 (MeV) 


1507 ± 5 


1507 


0.002 


Mf 3 (MeV) 


1714 ± 5 


1714 


~ 


J- f 2 ->TTTT (MeV J 


oo.U ± 4.0 


at) .84 


U.iOO 


1 h^KK (MeV) 


9.4 ± 1. ( 




n oqpi 
U.zoU 


r /2 ^ I)I( (MeV) 


5.6 ± 1.3 


4.89 


0.297 


1 ./WW 1 f 3 -,KK 


U.2U ± 0.06 


r\ i on 
0.129 


1.416 


^fs^nv/^fa—fKK 


0.48 ± 0.15 


0.227 


2.842 




0.88 ±0.23 


1.086 


0.863 


Tao — >7r?7 / /Tao — >7rr? 


0.35 ±0.16 


0.291 


0.134 


rif^Kx (MeV) 


273 ±51 


45.05 


20.365 


(r /3 ) 2P (MeV) 


140 ±10 


142.95 


0.087 


Xtot 






26.330 



Although the fit results given in Table 10 are acceptable, the smallness of the glueball-quarkonia mixing is clear 
contrast to other phenomenological studies Q-Q] and lattice result 
The mixing matrix reads: 



/ (A) ee |/ (1370)> \ / 0.82 -0.07 0.57 \ ( \N) = \nn) \ 

\f 2 ) = |/ (1500)> = -0.57 ~ 0.82 \G) = \gg) . (69) 
\ |/ 3 ) = |/o(1710)) J \ -0.06 0.99 0.04 / \\S) = \ss) J 

In this solution the state /o(1710) is very close to a pure gluonic configuration, which is traced to the small mixing 
parameter /. The states /o(1370) and /o(1500) are in turn dominated by the quarkonia components, but with strong 
mixing between nn and ss. The decay rates for /1 = / (1370) and fs, = /o(1710) as compared to the WA102 data are 
listed in Tables 11 and 12. 



Table 11. Decays of f x = / (1370). 



Quantity 


Exp (WA102) 


Theory 


F/i ->«/^/i-"" t 


0.46 ±0.19 


0.73 


r fx — > nn 1 r fx — > 7r tt 


0.16 ±0.07 


0.13 


(T fl ) 2p (MeV) 


"small" 


99.47 


Table 12. Decays of f 3 = /o(1710). 


Quantity 


Exp (WA102) 


Theory 


^ h^KKI^- /3->7r7r 


5.0 ±0.7 


7.78 


r/ 3 _ >nn f3^>Kn 


2.4 ±0.6 


1.76 


^/a— hjjj' /r/ 3 — unr 


< 0.18 


1.00 


(T /3 ) 2P (MeV) 


" dominant" 


142.95 



Finally, the two-photon decay ratios resulting from solution IV read: 

r /l ^ 27 : r /2 ^2 7 : r /3 _2 7 : r a o^ 27 = 1 : 0.273 : 0.008 : 0.382. (70) 
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E. Discussion of solutions III and IV 



The study of the previous subsections shows that a direct glueball decay with large violation of flavor symmetry 
is feasible. This feature is based on an interpretation of first lattice results, which corresponds to the limiting case 
of a direct glueball decay as contained in the parameters c d and c^). This interpretation is in contrast to the 
phenomenological study of |2( and to the comments given in |24j . Furthermore, such a large value for the flavor 
symmetry breaking parameter c 9 m is not in agreement with large N c arguments. In any case, considering that the 
interpretation of lattice results is not unique, from a phenomenological point of view it is interesting to analyze the 
situation, where a strong, flavor symmetry violating glueball decay is present. 

A sizable direct glueball decay is, as already discussed, the result of the phenomenology given in 8]. However, in 8] 
the direct glueball decay pattern is flavor blind, which would correspond to a large c d but to a suppressed value for 
c a m . In this sense the solutions III and IV differ from the analysis of |8j. 

Care should also be taken when considering the two-photon decay in this scheme: if the two-pseudoscalar amplitudes 
are sizable, the same can also be expected for the transitions into two vector mesons (although not studied here). 
Invoking vector meson dominance sizable corrections to the two-photon final state are expected. The use of the limit 
cf = in the present case is therefore questionable. 



V. CONCLUSIONS 



In this paper we analyzed the two-pseudoscalar and the two-photon decays of the scalar states between 1-2 GeV in 
the framework of a chiral Lagrangian, where the glueball has been included as a flavor-blind composite mesonic field 
with independent couplings to pseudoscalar fields. 

In a first step we have set the glueball-pseudoscalar couplings to zero and performed a fit to the accepted averages of 
PDG2004 0]. We find two possible solutions (I and II), which, apart from the underestimate of the Kq — > Kn mode, 
show good agreement with the data. The solutions I and II differ in the bare isoscalar masses, in the mixing matrix (in 
scheme I) the state /o(l 500) has the largest gluonic amount, while in scheme II the state /o(1710) has the main gluonic 
component), in some predictions concerning other decay modes, which have been compared to the experimental results 
of |33| . From a phenomenological point of view, there is no striking hint for a direct glueball-pseudoscalar coupling. 

We then enlarged our analysis by including the direct glueball decay parameter (c d ) in the flavor symmetry limit 
in the fit. A small value for this parameter is obtained, which in turn confirms the suppression of the direct glueball 
decay in agreement with large N c arguments and with 0, Q] , but contrary to the study of 8] . In a last step we also 
included the second glueball decay parameter (c^J in the fit, which involves flavor symmetry breaking in the direct 
glueball decay. The minima in the fit are less pronounced, therefore we utilized the lattice results of |2J| to determine 
the direct glueball decay parameters c d and c 9 m . The lattice data are interpreted such as they are matched by the 
direct glueball decay mechanism, although this procedure might be a limiting case. The resulting fits also generate 
a good description of the data, where either the /o(1500) or the /o(17IO) contain a dominant glueball component. 
These solutions however should at present taken with some care and require further input either from an enlarged, 
reliable database or lattice constraints. For this reason our preferred solutions are I and II presented in section III. 

Although the presence of a strong direct glueball decay cannot be verified directly, the presence of a sizable glueball- 
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quarkonia mixing is essential to be in accord with the data. The magnitude is different in the two proposed solutions 
(smaller i n I, largest in II), but is in line with other models and in magnitude qualitatively consistent with lattice 
results [iLll. 

The starting point of the Lagrangian has been outlined in J^JJ. In the present work we added the glueball degree 
of freedom, both for mixing and decays. As a result of the fit, we also considered deviations from the large N c limit. 
We find that large N c arguments are still useful as a guideline in the scalar sector as well. Also a small deviation from 
the GMO octet mass relation is found. 

A direct isoscalar nn-ss quarkonia mixing has been introduced in the theoretical analysis. Instanton solutions of the 
QCD vacuum are believed to generate strong flavor-mixing in both the isoscalar and the scalar sectors. The presence 
of such mixing is established in the pseudoscalar nonet, while it is still an open question in the scalar-isoscalar mesonic 
sector. Our two proposed solutions differ in this point: while in the first one this flavor mixing sensibly affects the 
results, in the second solution it turns out to be negligible. 

The problem of the Kq — ► Kit mode has been discussed; although no final statement to this puzzle can be said, we 
compared predictions of various approach and attempted to highlight the difficulty in the scalar sector. 

Although many experimental results can be reproduced, and the presence of a scalar glueball and its mixing 
with scalar quarkonia explains many features of the scalar meson spectroscopy, further work, both theoretically and 
experimentally, is needed to rule out some mixing scenarios in favour of others. 
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APPENDIX A: ASPECTS OF THE CHIRAL LAGRANGIAN 

We use the following phase conventions for the pseudoscalar P = {n^, tt°, K^jK , K°, r]°,rj 8 } and scalar 
S = {a^ , a®, Kq^, Kq°, Kq°, S°, S s } meson fields (neglecting the mixing of the third and the eighth component): 

7T ± = -^(PlT^2), 7r°=P 3 , K± = -^=(P 4 T»P 5 ), 
K° = -^=(P 6 - iP 7 ) , & = -±=(P 6 +iP 7 ), rj° = Pb, V 8 = Ps, 

clq = -y=(SiTiS2) , a,Q = S 3 , — 7^(^4 T 1S5) , 

K*° = ^(S 6 - zS 7 ) , = -L (5 6 + zS 7 ) . 

APPENDIX B: TWO-BODY s -> pip 2 AND s -> 77 TRANSITIONS (MATRIX ELEMENTS AND DECAY 

WIDTHS) 

1. Scalar-isoscalar strong decays 

The strong decay widths of the scalar states are derived at tree- level from the following term of the Lagrangian 

C-tZy = 4 (S u» un+c s m {S X+ ) (« M U ") + ^ G ( X+ ) • (Bl) 

The decay expression for the scalar-isoscalar states |z) with i — f%, fa, into mr, KK, r/rj and 7777' are given by the 
following expressions: 

V -V +r -*!r _ 3AV 2 (M 2 ,M 2 ,M 2 ) 

i i^TT-rr — 1 i^7r+7r- T - 1 - i— ►•7r°7r° — 2 i^7r+ir _ — 32 7T A'/? I j "i7r+ir _ I ! l,- DZ / 

r -r + r -,r _ AV 2 (M 2 ,M|,M|) 

M-siO? — 1 i-^-if+K- + 1 i^X»Jf» — 2 1 i^K+K Sir M 3 \ M iK+K-\ , \&<1) 



Ti—urn — 32ttM 3 l^wl ' (^4) 



z — ^7777' 



AV 2 (Af 2 ,M 2 ,M 2 ; 
32ttM? 

AV 2 (M 2 ,M 2 ,M 2 ) 



M lw ,| 2 (B5) 



16 7rMf 

where X(x,y, z) is the Kallen triangle function: 

A(ar, y, z) = x 2 + y 2 + z 2 - 2xy - 2yz - 2xz . (B6) 
The matrix elements M i7r +^~, M iK + K -, Mi m and Mi m i are given by 

= ~ 2M 2 ]c* d + 2M^} - |^|{[^ - 2M 'l^ + 2M 2 cf„}, (B7) 
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M^ 



K+K- 



2 5 



F 2 %/3 
2 c 



|[M? - 2M 2 K ]c 9 d + 2M 2 K clY 



[M? - 2M 2 ] |s iA r sin 2 5p + B 4S cos 2 <5 P n/2 j 



4 cf 



f> ^ _ | A/ 2 B lJV sin 2 S P + [2 Ml - M 2 }B lS cos 2 5 P >/2 
|^| {eg [A/ 2 - 2A/ 2 ] + 2 cf„ (a/ 2 + 2 [Ml - M 2 ] cos 2 * P ) } , 

S j^^ d [Mi - Ml - M 2 ][B iN - 5,5 V2] 

2c s m {Ml B lN + {Ml - 2M|) B iS V2] - c? n S iG [M| - Af 2 ]| , 



(B8) 



(B9) 



(BIO) 



where 5p = Op — 6 P and 9 P is the ideal mixing angle with sin#p = l/\/3; the quantities Bij are the elements of 
mixing matrix relating physical states i — f\,fa, /a and bare states j = N,G,S (see definitions in Eas. lTlj l and PTJll. 

2. Isovector and isodoublet strong decays 



The decay rates for ao(1450) into KK, iri] and 7tt/ are: 

1 



T o 



a°^K°K° 



X 1 / 2 M 2 ,M 2 ,M 2 , 



(Bll) 



O — »7T?J 



AV 2 (Af 2 o ,M 2 ,M 2 ) |2 



16ttM3 



(B12) 



A 1 / 2 (M 2 , M 2 ,M 2 ) 



16 7T M„ 3 



M 



The matrix elements M a + K+Ro , Af Q + 7r+j) and M a + ?r+r; , are given by: 



(B13) 



M„+ 



M„ 



M, 



--^([M o 2 -2M|]c^ + 2M|c r s 



([< ~ M 2 - MlY d + 2M 2 c 
_ cosW2 _ K _ + mX 



(B14) 
(B15) 
(B16) 
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The decay rates for Xq(1430) into Kir and Krj are (considering the isodoublet {Kq , Kq °}) 



2 ^K* + ^K°ir + 



32 7T Mj.. 



(B17) 



16ttML 

^ v n 



The matrix elements M K * + K07T + and M K *+ K+ are given by 



(B18) 



(B19) 



M K*+K + r, 



UMl, - M| - M 2 ] 4 + [3A4 - Ml] c s m 
+|^| ([M|. -Ml- Ml] 4 + [Ml + Ml] < 



(B20) 
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